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I. INTRODUCTION 

Numerical relativity simulations of binary black holes 
have matured dramatically in recent years. Starting with 
the first complete simulations of binary black hole merg- 
ers [IH3], a large number of papers on mergers of bina- 
ries with varying mass ratios and black hole spins have 
appeared. Some of the results of these simulations, in- 
cluding the surprisingly large recoil speed of the merger 
remnant for certain spin orientations (e.g. [H [5]), have 
also triggered numerous studies of the astrophysical con- 
sequences of these findings. 

Many numerical simulations of black holes adopt the 
Baumgarte-Shapiro-Shibata-Nakamura (BSSN) formula- 
tion of Einstein's equations [SJ [7] together with moving- 
puncture coordinates 0|3] (see also [S]). The latter con- 
sist of the 1+log slicing condition for the lapse [9j and 
the F-driver gauge condition for the shift 10J . The role 
of these coordinates in stabilizing the numerical simula- 
tions has been clarified by considering the late-time be- 
havior of spatial slices of the Schwarzschild spacetime 
when evolved with moving-puncture coordinates |11H13| . 
In particular, these studies showed that, at late times, 
these slices form a "trumpet" geometry, meaning that 
the slices asymptotically approach a finite areal radius 
and never reach the spacetime singularity (see Fig. 2 in 
[13) for an embedding diagram that motivates the name 
of this geometry). 

Several groups have also started to simulate black holes 
in higher dimensions (e.g. [T4lI22j ). Some of these stud- 
ies aim at exploring the rich geometric structure of black 
holes in higher dimensions (see, e.g., .23)), while others 
are motived by speculations that miniature black holes 
might be created in high-energy collisions in particle col- 
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liders (e.g. the LHC), a scenario that requires additional 
spacelike dimensions (see, e.g., [21] for a review). 

Motivated by the success of the BSSN formalism with 
moving-puncture coordinates in D = 4 spacetime dimen- 
sions, several of the above simulations adopt this method 
for D > 4 as well. In this paper we analyze the late- 
time behavior of spatial slices of spherically symmetric 
solutions in D > 4 spacetime dimensions when they 
are evolved with moving-puncture coordinates. More 
specifically, we construct stationary 1 + log slices of the 
Schwarzschild-Tangherlini spacetime [25] , the spherically 
symmetric vacuum black hole solution for D > 4. We 
demonstrate that, as for D = 4, these slices display a 
trumpet geometry. Parameterizing the 1 + log slicing 
condition with a parameter n, we show that we recover 
the results of [18] for maximal slicing in the limit n — > 00. 
We also recover the results of [13] for D = 4. Finally, we 
perform numerical simulations for D — 5 and demon- 
strate that these simulations indeed settle down to the 
trumpet slices. 

This paper is organized as follows. In Section |TT] we 
introduce the Schwarzschild-Tangherlini spacetime |25j 
and introduce a "height function" that we will use for 
coordinate transformations in the following sections. In 
Section |III| we construct a family of maximal slices 
of the Schwarzschild-Tangherlini spacetime, reproducing 
earlier results of |18j . and identify a special member of 
this family that displays a trumpet geometry. In Sec- 
tion |IV[ we derive the main result of this paper, a sta- 
tionary 1 + log trumpet solution for the Schwarzschild- 
Tangherlini spacetime. We recover the corresponding 
D = 4 results from [13] and the maximal slicing results 
of |18j in the appropriate limits. In Section[V]we perform 
dynamical simulations for D — 5 in spherical symmetry 
and demonstrate that, at late times, these simulations 
indeed settle down to the trumpet slices derived in the 
previous sections. We conclude with a brief summary in 
Section ED 

Throughout this paper we adopt units in which the 



2 



speed of light is unity, c = 1. However, we do not set 
to unity the gravitational constant G of a £>-dimensional 
spacetime, since keeping G in the analytical expressions 
makes their units more transparent. 



II. SCHWARZSCHILD-TANGHERLINI 

The generalization of the Schwarzschild spacetime, 
which describes spherically symmetric vacuum solutions 
in D — 4 spacetime dimensions, to higher dimensions D 
is the Schwarzschild- Tangherlini solution [55]. Adopting 
the notation of [23], this solution can be written as 

ds 2 = -f dt 2 + f^dR 2 + R 2 dn 2 D _ 2 , (1) 

where we have defined 



Eliminating ip we now find 



(2) 



&-D-2 — 



(3) 



(4) 



f<> ~ 1 R D-3 ' 

and where the mass parameter /i is given by 

16nGM 
M ~ (D - 2) n D _ 2 ' 

The area of a unit D — 2 sphere is 

2 7 r^- 1 )/ 2 

r((D-i)/2)' 

and the line element on this sphere is 

d ^D-2 = d xl + sin 2 X2dxl H 1" ( J[ sin 2 Xi^j dx D -v 

(5) 

where the xi (with £ = 2, 3, . . . , D — 1) are the angles on 
the D— 2 sphere. The angles X2 and xd— l are often called 
8 and 0. We also point out that GM has dimensions 
of (length) ~ 3 . For D = 4 we recover /x = 2GAf and 
f2 2 = 47r, as expected. 

In the line element (fTl) , R is the generalization of the 
Schwarzschild (or arealj radius. In numerical applica- 
tions, it is often convenient to express the metric in terms 
of an isotropic radius r, so that the spatial part of the 
metric can be transformed to cartesian coordinates very 
easily. To do so, we set the spatial part of the line element 
([!]) equal to its isotropic counterpart, 



dR 2 
1 - ii/R D - 



R 2 dn 2 D _ 2 = v/ /(I3 - 3) (d. 



r- + r 2 dQ 2 



D-2 i 



(6) 

where the exponent on the conformal factor ip has been 
chosen for convenience. From the identification Q we 
see that 



and 



i?2= ^4/(D-3) r 2 ; 

= v 4/(Z) - 3 W 2 . 



dR 



1 - ujR D 



(7) 



(8) 



± 



dr 
r 



R D ~ 3 dR 



y jR2{D-2) _ ' 

Integrating both sides of the equation we obtain [2^] 

2/(15-3) 



(9) 



or, solving for R, 

( a \2/(c-3) 

ii = K 1 + 4^) 

Inserting this into ^ we find the conformal factor 

/' 



tp = I 



4r 



D-3 • 



(10) 
(11) 

(12) 



As expected, these results reduce to the usual 
Schwarzschild expressions for D — 4. In Section [V] we 
will adopt the above expressions for D = 5 as initial 
data. 

In order to explore alternative slicings (or foliations) of 
the Schwarzschild- Tangherlini spacetime in Sections |III| 
and IIVI wc now introduce a new time coordinate 



t = t + h(R), 



(13) 



where h (R) is the height function [8l [27] . As we will see 
in the following Sections, different slicing conditions re- 
sult in different ordinary differential equations for h (R) . 
With the new time coordinate, the line element ([I]) be- 
comes 



2 h /2 



ds 2 



-fodtr + 2f h'dtdR - 



i - m 

fo 



-dR 2 + R 2 dn 



D-2 

(14) 

where h'(R) = dh/dR. 

We can write a general line element in (D — I) + f form 

as 



ds 2 



2 dt 2 + 7ij (dx l + P'dt) (dx 1 + fidt) , (15) 



where a, f3\ and jij are the lapse, shift vector, and spa- 
tial metric, respectively, and where the indices run 
over all D — I spatial coordinates. In spherical symme- 
try (which is preserved by the coordinate transformation 
(13 1) the only non- vanishing component of the shift is 
the radial component (3 R . Comparing terms in the line 
elements (14 1 and (151, we can identify the lapse as 



fo 



the shift as 



Wo 2 /*' 2 ' 

fph' 
l-/ 2 /*' 2 ' 



(16) 



(17) 
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and the spatial metric as 



(a- 2 
R 2 




7ij= 



R 2 sin 2 X2 



V 







R 2 l\Z~2 2 sin 2 X// 



The determinant 7 of the spatial metric is 

D-2 



7 = a 



J i? 2 (^ 2 > J] (sin 2 Xi 



D-l-l 



1 = 2 



We define the extrinsic curvature so that 



dtlij 



-2aKij + Difa + I), 1, 



(18) 
(19) 

(20) 



where Di is the covariant derivative operator associated 
with the spatial metric. For time-independent spatial 



metrics like (18 1, the left-hand side of (20 1 vanishes and 
we obtain 



2^- = ^(A& + £>iA). 



(21) 



Taking the trace of this equation we find that the mean 
curvature is given by 



K = 



1 



V2 dR 



cry 



(22) 



This equation is our starting point for imposing maximal 
slicing in Section ITTT] and 1 + log slicing in Section [TV] 



III. MAXIMAL SLICING 

Maximal slicing is defined by requiring that the mean 
curvature vanish, 



K = 0. 



(23) 



Results for maximal slices of the Schwarzschild- 
Tangherlini spacetime have already been presented in 
[15] ; the details included here are for the sake of com- 
pleteness and reference in later sections. In Section [lII A[ 
we derive a family of time-independent maximal slices 
of the Schwarzschild-Tangherlini spacetime for a general 
number of spacetime dimensions D > 4, and we special- 
ize to D = 4 and D = 5 in Sections |mB| and pTC] 
respectively. 



General treatment 



For maximal slicing, equation (22) reduces to 
d fR D - 2 p Rs 



dR 



= 0. 



(24) 



Eliminating a and f3 R with the help of equations ( 16 ) 
and (17) we obtain a first integral 



R 



D-2 



fa 



1 - f 2 h' 2 



1/2 



foti = C, 



(25) 



where C is a constant of integration. It is convenient to 
write ( 25 1 as 



f h' 2 = 



C 2 

f R 2( ~ D - 2 ) + C 2 ' 



(26) 



which we can then substitute back into equations ( 16 ) 



and (17) to find the lapse 

a = f(R;C) 



the shift 



R Cf (R; C) 



R D ~ 2 



and the spatial line element 

dl 2 = r 2 (R; C) dR 2 + R 2 dn 2 D _ 
Here the function f(R; C) is given by 



f(R;C) 



1 



c 2 



R D-3 R 2(D-2) 



1/2 



(27) 



(28) 



(29) 



(30) 



The one-parameter family of spherically-symmetric, 
time-independent maximal slices of the Schwarzschild- 
Tangherlini spacetime is now parameterized by the con- 
stant C . For sufficiently small C the slices end at a radius 



R at which the lapse vanishes; from equation (27), this 



location is given by the largest root (see [TB]) of the equa- 
tion 



R, 



2{D-2) 




D-l 



c 2 = 0. 



(31) 



Two particular members of this family deserve spe- 
cial mention. For C = the height function h must 
be constant, the spacetime is therefore sliced by slices of 
constant Schwarschild-Tangherlini time t, and we recover 
the metric ([!]). The other member of the family that we 
will be interested in is that for which the slice ends with 



a double-root of the squared lapse (27). As we will show 
below, this choice singles out a trumpet slice. In Sec- 
tion [V] we will see that these slices act as "attractors" in 
dynamical moving-puncture simulations. 

Setting both a 2 and the first derivative d(a 2 )/dR to 
zero, we find that this double root occurs at the radius 



Ro 



H{D-1) 



2 (D-2) 
for the value of C given by 



C 2 



D-3\ ///(D-l) 



D-l J \2(D-2) 



l/(D-3) 



2(D-2)/(D-3) 



(32) 



(33) 
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where the tilde denotes the special value for the trumpet 
slice (see also [18]). 

Before closing this section we analyze the asymptotic 
properties of these slices in a neighborhood of Rq. Given 
the spatial metric (181, we can compute the proper dis- 



(21 1, respectively. Also, for the special value of C, which 
3\/3(GM) 2 



now becomes 



C 



(38) 



tance A between a point on the limit surface, (Ro,Xi) 
and a point at (Ri,Xi) from the integral 



the transformation to isotropic coordinates can be carried 



out by integrating equation ( 36 1 (see [30] 



A = 



fll dR 
Ro a ' 



(34) 



By definition, the lapse a vanishes at Ro. Whether or not 
this integral is finite therefore depends on the behavior 
of a in the neighborhood of Ro. From its definition in 
Eqs. (27 1 and (30), we see that the square of the lapse 



can be expanded as 



C. Five-dimensional spacetimes 

For five-dimensional spacetimes, D = 5, we recover 
the results of [18J. We note, however, a difference in 
notation: our R is their r and vice versa. With fl 3 = 2tt 2 
and /i = 8GM/(3ir) we now find the lapse 



A 1 {R-Ro)+A 2 {R-R f 



(35) 



a = / (R; C) = 1 



8GM 



c 2 \ 1/2 



3irR 2 R e j 



(39) 



For generic values of C, R = Ro is a single root of a 2 
and hence Ai is nonzero. Therefore the lapse beh aves 
like a ~ (i? — i?o) 1//2 near Rq and the integral in Eq. (34) 



With the lapse known, the shift, spatial metric, and ex- 
trinsic curvature can again be calculated from (28), (18), 



is finite. Thus, for generic C, the proper distance to the 
limiting surface at Ro is finite. 

For the special slice C = C, on the other hand, R = Ro 
is a double root of a 2 and therefore A\ vanishes. The 
lapse then behaves like a ~ (R — Rq) near Rq and the 



and (21), respectively. The trumpet geometry is realized 
for 



- _ 2"/2 /GM\ 3/2 



integral in Eq. pih diverges. The proper distance to the and its limiting surface is at 



limiting surface at Rq is infinite, even though the areal 
radius is non-zero and finite (see |18j for an alternative 
derivation of this property). An embedding diagram of 
the slice would result in a figure similar to that in Fig. 2 
of [13] ; given the appearance in an embedding diagram 
these slices are referred to as "trumpet" slices. 

As in Section [IT] it would be useful to transform these 
slices to isotropic coordinates. Instead of equation 
we now have 



Ro 



4 ja\j 



(40) 



(41) 



For the trumpet slice, equation ( 36 1 can again be inte- 
grated analytically (see 
isotropic coordinates. 



]) to obtain the solution in 



± 



dr 



1 dR 



R D ~ 3 dR 

y/R^D-2) _ ^ R D- 



C 2 

(36) 

For C = C, the expression under the square root on 
the right hand side has a double root at R — Rq, which 
simplifies the integral for both D = 4 and D = 5. 



B. Four-dimensional spacetimes 

For four-dimensional spacetimes, D = 4, we re- 
cover the well-known family of time-independent max- 
imal slices I27H55] . In particular, with f2 2 = 4-zr 
and fx — 2GM, the lapse becomes 



IV. STATIONARY 1 + LOG SLICING 



We now turn to stationary 1+log slicing 
(3 l diOt — naK, 



(42) 



where n is a constant. The most common choice for n is 
n = 2, but, following [T3], we will pursue a more general 
treatment. In particular, we will recover the maximal 
slicing results of Section III in the limit n — > oo. Mirror- 
ing Section |III[ we derive a family of time-independent 
1 + log slices of the Schwarzschild-Tangherlini spacetime 
for arbitrary D > 4 in Section IV A[ and specialize to 
D = 4 and D = 5 in Sections ITVl¥ a7d llVCl 



A. General treatment 



f(R;C)= i 



2GM 
R 



91 
i? 4 



1/2 



(37) 



Using equations (14), (15), and (22), we can rewrite 
(|42l as 



With the lapse known, the shift, spatial metric, and ex- 
trinsic curvature can be calculated from (28), (18), and 



— Oi = 

dR 7 1 /2 / 3fl dR 



V2fl« dR V P J 



(43) 
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It is convenient to use (16) and (17 1 to find 
(3 R = a^a 2 - f , 



(44) 



and to eliminate f3 R in ( 43 ) . This substitution results in 
'/ , , ' ■ " ' (45) 



dR a ~ U dR 



In (R^V^fo), 
which can be integrated immediately to yield 
a = const + n In [r (d ^ 2) a/ck 2 - f ) , 



or equivalently, 



/o + 



C 2 (n) e 2a ' n 

i?2(D-2) ■ 



(46) 



(47) 



Here C is again a constant of integration, but we point 
out that it now depends on the constant n. Equation 



(47) is a transcendental equation for the lapse a, and so- 



lutions, when they exist, can be found numerically. With 
the solution for the lapse in hand, the shift, spatial met- 
ric, and extrinsic curvature can be found from (44), (18), 



and ([21 >. 

As for maximal slicing, some values of C 2 (n) are 
special. Setting C 2 (n) = corresponds to the usual 
t = const slices, recovering the metric 0. Another so- 
lution of interest is obtained by requiring equation ( 45 ) 



to be regular for all a > 0. Following the method given 
in [13] for D = 4, we can use the regularity condition 
to determine the corresponding values C 2 (n) for D > 4. 
We begin by rewriting ( 45 1 as 



dR 



a 



(P-2)j?- 1 (o 2 -/ )-i/-) 
fa + not — a 2 



(48) 



If we require regularity when a > 0, the numerator and 
denominator must vanish at the same value of R. The 
numerator is zero when 



Q 



7o + 



2 (D - 2) ' 



(49) 



Substituting this value for the lapse and /o from ^ into 
the right hand side of ( 48 1 , we see that its denominator 
vanishes when 



R 



2 (D - 2) 



QP-3)V 
4n 2 (13 - 2) 2 



0. (50) 



This equation is quadratic in i? 15 3 , and so it is easily 
solved for the critical value R r for which the numera- 



tor and denominator of ( 48 ) vanish simultaneously. The 
positive real root is 

Rr = 



(D-l)u + uJ4(D-3) 



(D-iy 



4(D-2) 



(51) 



The lapse at R c can be calculated from ( 49 ) to be 



4 (D - 3) 2 + (D - if n 2 -(D- 1) n 



\ yj 4 (D - 3) 2 + (D - if n 2 + (D- 1) n 



(52) 



Inserting equations (51) and (52) into (47) we now find 
C 2 (n) , defined as the special value of C 2 (n) that makes 
(48) regular, 



C 2 (n) = 



D-2 



D-3 



x \ D 



4{D-3f + (D-lfn 2 



D-l 
D-3 



n 



(53) 

For most of the work that follows we consider only the 
special slices with C = C. The inner boundary of these 
slices is defined by the location at which the lapse goes to 
zero. We will refer to this location as the throat, and will 
call its radius Rr,. To identify Rq, we first observe that, 



at least for sufficiently large R, equation (47) admits two 



solutions for the lapse a. For R — > oo, one of the branches 
of solutions approaches a = 1, while the other approaches 
a = — 1. This behavior is displayed in Fig. [T] where we 
graph a versus R for D = 5 and n = 2, for the special 
value C = C. We are interested in the "positive" branch 
of solutions, which approaches a = 1 asymptotically. As 
can be seen from Fig. [T] this positive branch has a root 
(with a = 0) at a smaller value of R than the negative 
branch. In general, then, we define the throat as the 
location Rq at which the positive branch of solutions for 
a vanishes. For D — 4 and D = 5, this corresponds to 
the smallest root of equation (47) with a = 0, 



R 



2{D-2) 




uRV- 1 + C z (n) = 



(54) 



(see also the discussion in [13 ). In the limit n — > oo, 
the critical value of the lapse a c — > and the slice ends 
in a double root, so this is still consistent with choosing 
the largest root in the case of maximal slicing. We will 
discuss explicit solutions for Rq (n) for D — 4 and D = 5 
in Sections IVB] and HVC] below. 

We can demonstrate that the special slice with C = C 
is a trumpet slice. To do so, we assume that close to Rq, 
a can be expanded as a power series 



R 



a 



5 

i? 



(55) 



where R = Rq + 5. Since a = for S — we must have 
m > 0. Inserting this into equation (47), with C = C, 



and expanding terms to first leading order, we obtain 

2m ^,2 < 




C 2 



n2(D-2) 



R 







^ nD-3 



6 

Ro 
2(2- D)- 



C 2 



2{D-2) 



R 



(I 



_5_ 

Rq 



(56) 
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FIG. 1: The two branches of the lapse a as a function of the 
radius R for D — 5, n = 2, and C = C. We are interested in 
the "positive" branch that approaches a = 1 at infinity. This 
branch has a root (with a = 0) at a smaller value of R than 
the "negative" branch. 



Matching exponents of S we identify m = 1, so that 
a ^ 5 = R — Rq near R = Rq. Therefore the proper 
distance from any point out sid e the throat to the throat 
at Rq, given by the integral ( 34 ), diverges. As before, this 



demonstrates that this slice features a trumpet geometry. 

In the limit n — > oo, the 1 + log slicing condition (42) 
approaches the maximal slicing condition K = 0, so we 
should recover the maximal slicing results from Section 
TTll From (|47| we have 



o 



/o 



R2(D-2) ■ 



consistent with (27 1. Similarly, taking the limit n 



in equation ( 53 1 we recover ( 33 ) . Finally, letting n 



in ( 51 ) results in 



R r 



mQp-i) 

2(13-2) 



l/(£>-3) 



(57) 
-> oo 

00 



(58) 



Substitution into ( 54 ) shows that this coincides with the 
location of the throat, so that R c = Rq in this limit, 
meaning that we have recovered Rq in the maximal slic- 
ing limit ( 32 ) . 

As in Section [TTJ it would useful to transform these 
slices to isotropic coordinates. While it is impossible to 
carry out the necessary integrations analytically, given 



the transcendental nature of equation ( 47 1 , the equations 



can be integrated numerically using the techniques pre- 
sented in IT51. 



FIG. 2: The throat location Rq as a function of the coefficient 



n in the 1 + log slicing condition ( 42 1 for D — 4 (dashed) 
and D — 5 (solid). The points are the corresponding values 
found by our numerical simulations described in Section [V] 
Specifically, we extrapolated the lapse a linearly to zero to 
find the throat location _Ro for any grid resolution, and then 
used Richardson extrapolation using our two highest-order 
resolution results to obtain the values shown here. 



B. Four-dimensional spacetimes 

For four-dimensional spacetimes, the above results re- 
duce to those of [15] . In particular, for D = 4, equation 
( 54 1 becomes a quartic equation 



Rl - fiRl + C 2 (n) = 0. 



(59) 



The real root consistent with the positive branch of the 
lapse is 



Ro(n) = » 1 + J 1 + Z-J 1 -Z+ - 1 
°^ 2 I 2 V 4 y 2 4^1/4 + Z 



where 



Z = 



4(l) 1/3 g 2 (n) 
F/it 4 



Y 



21/332/3 ' 



and 
Y 



9 + t/81 - 768C 2 (n) ^ 4 (n 



1/3 



(62) 

We graph these solutions in Fig. ^ (compare with Fig. 3 
of |13] , where these solutions are given only graphically) . 
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For comparison with [T3] , we note that ( 48 1 can be 
written as 



dR a R{R D - 3 (f + na-a 2 )) 



{-(D-2)R^f --R^ 2 f a 
■(D-2) R°- 3 a 2 ), 



D-2 el (63) 



which, for D = 4, reduces to 



d n (3GM -2R+ 2Ra 2 ) 

dR a ~ R(R — 2GM + nRa - Ra 2 ) ' 



(64) 



From this equation, the remaining results of Section 
II D 1 of H3| can be derived. 



C. Five-dimensional spacetimes 



For D = 5, equation (54 1 becomes a cubic equation for 
the Rl(n), 



Rl - fj,R^ + C 2 (n) = 0, 



(65) 



which can again be solved exactly. The root that is real, 
positive, and consistent with the positive branch of the 
lapse is 



Rq 



Jl/2\ 



1 / 27C 2 (n)\ 4tt 

2cos -arccos 1-^^ + y 



(66) 



1 



We show a graph of Rq in Fig- [2] 

With n 3 = 2tt 2 , pi = 8GM/(3tt), and f 
8GM/(3irR 2 ), the solution ([47} of the stationary 1+log 
slicing condition becomes 



of = 1- 



8GM C 2 (n) e 2a / r 



3ttR 2 



R 6 



(67) 



This determines a as a function of R. When solutions to 
this equation exist, they can be found numerically with 
a root-finding routine by starting with an appropriate 
guess on the positive branch of the lapse. With the lapse 
known, the shift, spatial metric, and extrinsic curvature 



can be calculated from (44), (18), and (21), respectively 



For some purposes, it is sufficient to find R (a) instead of 
a (R). The real, positive solutions of relevance here are 
closely related to the solution (66) for Rq (n): 



Ri {a) 



3(1 -a 2 ) 



, 1 / 27A\ 4tt\ 
1 2 cos I - arccos II — 1 + — ) + 1, 



(68) 



and 



R2 (a) 



3(1 -a 2 ) 



x 4 / 2 cos [ - ai'cci >s 



(69) 



where 



A = — f ) C 2 {n)e 2a ' n . 



(70) 



With C = C, the positive branch of the lapse is given 
by R\ (a) for < a < a c and R2 (a) for a c < a < 1. 
Solutions for D = 4 can be constructed similarly from 
solutions to equation |59| ). 

The critical radius (151 D is now 



2 (2GM 

Rc = 5 I r, 

6 \ wn 



1/2 



and the lapse ( 52 ) at R = R c is 



vl + n 2 — n 
Vl + n 2 + n 



(71) 



(72) 



Finally, we find C from equation ( 53 ) , 



CHn) = ^^(l + 1 -V^) 2 e- 2 ^ (73) 

As before we can recover the maximal slicing results of 
Section III C by letting n — > 00. From equations ( 58 ) and 



( 66 ) we then obtain 



R c = Ro 



GM 



(74) 



which agrees with equation (41). 

Most numerical relativity simulations adopt n = 2. In 
this case equation (51) yields 



R c = ^\ 2 + ^ VGM sa 0.774132%/GM, (75) 
3 V K 



and from ( 53 1 we have 

128(2 + V5)W(^- 2 )/(^+ 2 ) 



C 2 (2) 



729tt 3 
0.0802483 (GM) 3 . 



The throat ( 66 ) is located at 



« 0.725474/i5 « 0.668392y / GM. 



(GM) 3 
(76) 

(77) 



For D = 5, the event horizon in the Schwarzschild- 
Tangherlini spacetime is located at -Reh = 
a/8GAJ/(37t) = 0.921318VGM. At the horizon, 
the lapse is 



ctEH ~ 0.454702. 



(78) 
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V. NUMERICAL WORK 

In this section we describe dynamical numerical sim- 
ulations of the Schwarzschild-Tangherlini spacetime for 
D = 5. More specifically, we adopt as initial data the 
t = const slices of the Schwarzschild-Tangherlini space- 
time in isotropic coordinates, as described in Section [TT} 
and evolve these data with moving-puncture gauge con- 
ditions. We will see that the trumpet solutions derived 
in the previous sections indeed act as "attractors" in 
dynamical simulations, meaning that dynamical simula- 
tions settle down to these solutions at asymptotically late 
times, even when they start with very different initial 
data. 



which the conformal factor is given by equation (12 1. If 



these data were evolved with the Killing lapse and shift 
that can be identified from the metric (JTJ), the metric 
would remain time-independent. Instead, we set 



a 



1. 



(80) 



initially (at t = 0), and subsequently evolve the lapse 
and shift using moving puncture coordinate conditions 
O [3], namely the 1+log slicing condition for the lapse 
[H] and the T-driver condition (53]. We will consider two 
different flavors of the 1+log slicing condition, namely an 
"advective" and a "non-advective" version. 



1. 1+log slicing 



A. Numerical Method 

Our numerical code is based on the third-order finite- 
difference code of [31]. This code implements the BSSN 
equations [5] [7] for D = 4 in spherical symmetry, and 
imposes spherical symmetry with the help of the "car- 
toon" method [32] . Details of the implementation of this 
method can be found in Appendix A of [31j . Several mod- 
ifications had to be implemented to adopt this code to 
five spacetime dimensions. Clearly, all indices of spatial 
tensors had to be extended to four spatial dimensions, 
and the interpolation for the cartoon method had to be 
adjusted for the extra dimension. In addition, the finite 
differencing stencil was updated to account for the fourth 
spatial dimension. Furthermore, even though Einstein's 
equations take the same form for all spacetime dimen- 
sions D, the BSSN equations for D — 5 are different from 
the corresponding equations in D — 4 (see, e.g. [TT])- The 
reason for this is the use of tracefree tensors in the BSSN 
formalism. In particular, consider the decomposition of 
the extrinsic curvature into its trace K and its trace- 
less part A 



The "advective" 1+log slicing condition is given by 



131 



Ki. 



D 



(79) 



Since the BSSN equations are formulated in terms of A%j 
and K instead of Kij , factors that depend on D appear 
in several places. Following [T7], we also use the vari- 
able x = V' -2 instead of the conformal factor ip in our 
numerical simulations. 

We tested features of our code that are specific to D = 
5 with the help of a number of test problems, including a 
spherically symmetric standing wave. We also monitored 
constraint violations during our simulations and verified 
that they converged to zero at the expected rate. 



B. Numerical Results 

In all our dynamical simulations we start with data 
on a slice of constant Schwarzschild-Tangherlini time, for 



(d t - (3 l di)a = -naK, 



(81) 



where n is again a constant. At late times, when the so- 
lution settles down and becomes time-independent, this 
condition reduces to the stationary 1+log condition (42 1. 



Our numerical results demonstrate that simulations 



with the advective 1+log slicing condition (81) indeed 



settle down to the stationary 1+log results of Section 
IV C In the following we show results that were ob- 
tained with = 5001 uniform gridpoints (not count- 
ing the buffer points), with the outer boundary imposed 
at r max = 50%/ GM. We also chose n = 2 for all re- 
sults shown in this section. (We performed simulations 
for other values of n and N to obtain the results for the 
radius of the throat Rq shown in Fig. [2]) 

In the following figures we graph several quantities in 
terms of the areal radius R, which simplifies comparison 
with the analytical results for the late-time asymptotic 
solutions. In terms of the variables evolved in our code, 
we compute R from R = y^eoJXi where jgg is the 99 
component of the conformally related metric and x = 

The initial data for our dynamical simulations are 
given in terms of an isotropic spatial metric. However, 
this spatial metric is evolved in time, and does not re- 
main isotropic (see Figs. 35 and 36 in [UJ), so that we 
cannot compare the conformal factor in the code directly 
with the conformal factor that we would obtain by trans- 
forming the analytical results of the previous sections 
into isotropic coordinates. Instead, we compare analyti- 
cal and numerical values for the metric component 7^ 
in Fig. [3] We compute numerical values for 7^7? from 



7RR 



X (dR/dr) 



(82) 



where the rr component of the conformally related 

metric and the derivative dR/dr is calculated numerically 
from the areal radius R given above. 

We show the evolution of the numerical 7^ at sev- 
eral different times, starting with the initial data ( 12 ) 
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FIG. 3: The metric component 7_r_r as a function of the 
areal radius R for a dynamical evolution of the D = 
5 Schwarzschild-Tangherlini spacetime with the advective 
1+log slicing condition (81 1 for n = 2. We show ^rr at a 



FIG. 5: Same as Fig.^ but for the shift f3 R . 



number of different times, starting with the initial data ( 12 1 
t — 0, and ending at time t — 42V GM. At late times, our 
numerical results for ~/rr agree very well with the analytical 
results of Section IIV CI The latter are included as the solid 
line, which overlaps with the line representing the numerical 
results at t = 42VGM. 
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FIG. 4: Same as Fig. [3] but for the lapse a. 



at t = and ending at a time t — 42V GM before the 
results shown in the figure are affected by the presence 
of the outer boundary. (In order to avoid double-valued 
functions, we show the numerical initial data only for 
isotropic radii outside the event horizon.) Also included 
are the analytical results obtained in Section IV C We 
see that "/rr evolves for a certain period of time, and then 
settles down to the expected time-independent solution. 

In Figs. [4] and [5] we show similar results for the lapse 
and the shift. The numerical values of the areal shift /3 R 
are calculated using (3 R = (dR/dr)/3 r . Both the lapse and 
shift start with their initial values ( p0| , go through a dy- 
namical phase, and then settle down to their asymptotic 
values. At these late times we find excellent agreement 
with the stationary solutions that we computed analyti- 
cally in Section |"lV C| 



2. Maximal slicing 

We also performed numerical simulations for the "non- 
advective" 1+log slicing condition, for which the advec- 
tive term in (81 ) is dropped, 



d t a 



-naK. 



(83) 



Time-independent solutions must satisfy the maximal 
slicing condition, K — 0. If dynamical simulations with 
this slicing condition settle down to a time-independent 
solution, then the late-time asymptotic solution must be 
given by a member of the family of maximal slices that we 
discussed in Section [nil As it turns out, dynamical simu- 
lations settle down to the "special" member with C = C, 
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FIG. 6: The metric component 7_r_r as a function of the 
areal radius R for a dynamical evolution of the D — 5 
Schwarzschild-Tangherlini spacetime with the non-advective 
1+log slicing condition (831 for n = 2. We show 7_r_r at a 



number of different times, starting with the initial data ( 12 1 
0, and ending at time t 



t — 0, and ending at time t — 42v GA/. At late times, jrh 
approaches the trumpet member of the family of maximal 
slices discussed in Section Till CI The latter is included as the 
solid line, which overlaps the line representing the numerical 
results at t = 42VGM. 



VI. SUMMARY 

We study maximal and stationary 1+log slices of the 
Schwarzschild-Tangherlini spacetime for D > 4 space- 
time dimensions. We use a height function to introduce 
a coordinate transformation away from slices of constant 
Schwarzschild-Tangherlini time, and find families of both 
maximal slices and stationary 1+log slices. 

For the maximal slices, which were previously consid- 
ered by [IS], we identify one "special" member that dis- 
plays a trumpet geometry. For the stationary 1+log slices 
we impose a regularity condition following the methods 
of [13 . This regularity condition singles out one partic- 
ular member, which also displays a trumpet geometry. 
We allow for a free parameter n in our 1+log slicing con- 
dition, where the limit n —¥ oo corresponds to maximal 
slicing. We demonstrate that our results for 1+log slicing 
reduce to the maximal slicing results of |18j in this limit, 
and we also show that we recover the results of [13] for 
D = 4. 

Finally, we perform numerical simulations for spher- 
ically symmetric black holes in D = 5 spacetime di- 
mensions. We start with data on a slice of constant 
Schwarzschild-Tangherlini time, and evolve these data 
with moving-puncture gauge conditions. Our results 
demonstrate that, as for D = 4 spacetime dimensions, 
the dynamical simulations settle down to the trumpet 
slices, which can be regarded as "attractors" for moving- 
puncture simulations. 
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